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ABSTRACT 

By simulating jet-inflated bubbles in cooling flows with the PLUTO hydro- 
dynamic code we show that mixing of high entropy shocked jet's material with 
the intra-cluster medium (ICM) is the major heating process perpendicular to 
the jets' axis. Heating by the forward shock is not significant. The mixing is 
very efficient in heating the ICM in all directions, to distances of ~ 10 kpc and 
more. Although the jets are active for a time period of only 20 Myr, the mix- 
ing and heating near the equatorial plane, as well as along the symmetry axis, 
continues to counter radiative cooling for times of > 10 8 yr after the jets have 
ceased to exist. We discuss some possible implications of the results, (i) The 
vigorous mixing is expected to entangle magnetic field lines, hence to suppress 
any global heat conduction in the ICM near the center, (ii) The vigorous mixing 
forms multi-phase ICM in the inner cluster regions, where the coolest parcels of 
gas will eventually cool first, flow inward, and feed the active galactic nucleus to 
set the next jet-activity episode. This further supports the cold feedback mech- 
anism. (Hi) In cases where the medium outside the region of r ~ 10 kpc is not 
as dense as in groups and clusters of galaxies, like during the process of galaxy 
formation, the forward shock and the high pressure of the shocked jets' material 
might expel gas from the system. 



1. INTRODUCTION 



The intra-cluster medium (ICM) in cooling flow (CF) clusters and groups of galax- 



ies is heated by a negative feedback mechanism (e.g., iBinnev fc Tabor! Il995l : iFarage et al.l 



Birzan et al. 2011 



e.g., 


Done: et al. 


Gitti et al. 


2012 



20 id b 'Sullivan et~aJ l201l[ baspari et aD l2012al lbl: 

for recent papers and references therein). Examples 
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of bu bbles include Abel l 2052 (IBlanto n et al 



5044 (IDavid et a 



( iMachacek et al 



2009), HCG 62 (iGitti et al. 



201lh NOG 5813 (IRanda 



2011 



Abell 4059 (IHeinz et al.ll2002h . NGC 4 636 (IBaldi et al 



2009 iDavid et al-KOllt ). and RBS 797 flSchindler et al. 



1 et al. 



), NGC 633 8 (IPandge et al 



201 0h, Hyd ra A (TWise et al.ll2007l). NGC 584 6 



201 



2012b, NGC 



2001 



), A 2597 (IMcNamara et all 12 001) . 
20091). NGC 5044 (lOastaldello et al. 



Cavagnolo et aU2011uDoria et al 



2012) 



Wide bubbles very close to the origin of the jets (the AGN), e.g., as in Abell 2052, that 
are termed 'fat bubbles', can be inflated by jets that do not penetrate through the ICM. 
Instead, they deposit their energy relatively close to their origin and inflate the fat bubbles. 
Slow massive wide (SMW) jets can inflate the fat bubbl es that are observed in many CFs, 
in clusters, groups of galaxies, and in elliptical galaxies ([Sternberg et al. 1 120071 ). The same 
basic physics that prevents wide jets from penetrating through the ICM holds for precessing 



jets ([Sternberg fc Soker 



jets t o the medium (IBriiggen et al. 



2008a 



Fa 



ceta-Goncalyes et al. 201ol). or a relat i ve motion of the 



2007 



Soker 



2009: 



Morsony et al.ll2010t iMendygral et al. 



20121 ). In the present study we will inflate bubbles by SMW jets, but our results hold for 
bubbles inflated by precessing jets or a relative motion of the ICM as well. If the jets 



elongated and/or detached from the center bubbles are formed (e.g., 


Basson & Alexander 


2003; 


Omma et al. 


200' 


1; 


Heinz et al. 


2006; 1 


\/ernaleo & Reynolds! 


20061: 


A 


ouani Bibi et al. 


2007; 


Sternberg et al. \ 


2007; 


O'Neill & Jones 


2010; 


Mendvgral et al. 


2011, 


2012 


)• 



Vortices inside the bubbles and in their surroundings play major roles in the forma- 
tion of bubbles, their evolut ion, and their interaction with the ICM (e.g. IHeinz fc Churazov 
20051 ). lOmma et al.l (120041 ) find that a turbulent vortex trails each cavity , and that this 



vorte x contains a significan t quanti ty of entrained and uplifted material (also iRoediger et al. 



20071 ). ISternberg &; Sokerl (l2008bl ) find in their 2.5D numerical simulations that vortices 
inside bubbles can stabilize them against the Ray leigh- Taylor (RT) instability and can sup- 
press the Kelvin-Helmholtz (KH) instability on the surface of each bubble. Jet-excited 
shock s that interact with older bubbles can excite vortices that dissipate energy to the 



ICM ([Friedman et al.ll2012l ). The vortices cause semi-periodic changes in the bubble prop- 



(Sternberg & Sokerl \ 


2009). 


(Refaelovich & Soker 


2012) 



in the interaction of the jets and bubbles with the ICM. We concentrate on gas near the 
equatorial plane, i.e., ICM gas that does not reside along the jets' expansion trajectory. 

The heating of the gas perpendicular to the jets' axis need not be 100% efficient, as 
observations show that heating does not completely offset cooling (e.g., IWise et al.l I2004J ; 
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McNamara et all 12004; Oarke et all 12004 : iHicks & Mushotzkvl 12005 



Salome et al. and a moderate CF exists (ISokerll2004l ). Moderate 



Brcgma n et al.l 12006 



implies here that the mass cooling rate to low temperatures is much lower than the cooling 
rate expected without heating, but it is much larger than the accretion rate onto the super- 
massive bla ck hole (SMBH) at the center of the clu ster. The cooling mass is eit her forming 



stars (e.g., lO'Dea et al.l l2008t iRafferty et al.l 120081 ) . forming cold cloud s fe.g.. Edge et al. 



2010), accreted by the SMBH to maintain the cold feedback mechanism (jPizzolato fc Soker 



2010 ) . or is expelled back to the ICM and heated, when it is shocked or mixed with the hot 



jets' material. 

In the present paper we study the heating and expelling process of gas residing near 
the equatorial plane. By running cylindrically symmetric hydrodynamical simulations with 
the PLUTO code (section [2]) we study the flow structure (section [3]). We then examine the 
degree of mixing and dredge-up of gas from the equatorial plane vicinity (section HJ), and the 
heating processes (section E]). Our short summary is in section [HI 



2. NUMERICAL SET UP 



The simulations were performed using the PLUTO code (IMignone et al.ll2007l ). We use 
spherical coordinates in 3D, but with an imposed azimuthal symmetry around the 9 = (z) 
axis. Namely, a 2.5D grid where only the dependence on the (r, 9) coordinates is calculated. 
We will present the results in a plane <fi =constant, the meridional plane, which we take to 
be the (w, z) plane, where vo is in the equatorial plane and z is along the symmetry axis. 
The grid radial and azimuthal domains are r = 0.5 — 215 kpc, and 9 = — 90°, respectively. 
We focus on the relevant inner domain of r < 40 kpc; the large outer radius ensures that 
there are no boundary effects, at a low computational cost thanks to a radially stretching 
grid. 

On the r = 0.5 kpc inner boundary we enforce a jet outflow between angles — 70° 
for a limited time, and reflective boundary conditions in the angular zone 70° — 90°. In 
the simulations presented here the two jets power is P 2 j = 2 x 10 44 erg s _1 , their initial 
velocity is Vj = 9600 km s" 1 , and mass loss rate into the two jets is M 2 j = 7M yr _1 . Such 



mass i ve wide sub-relatiy i stic outflows are sup ported by recent observations (e.g.. iMoe et al. 



2009 



Dunn et al. 



2010 



Tombesi et al.l 120121 ) . At t = 20 Myr we turn the jet off, and 



apply reflective boundary conditions on the entire inner sphere. We set reflective boundary 
conditions on the boundaries 9 = [z axis) and 9 = 90° (the equatorial plane). The 
latter boundary condition mimics a symmetric opposite jet. The initial density profile has a 
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spherical symmetry, with a profile of the form (e.g., IVernaleo fc Reynolds! 120061 ) 



PM = Pc 



r 



(1) 



We take here r = 100 kpc, and p c 



10 25 g cm 3 . 



We start the simulation with an isothermal sphere at temperature T = 4 x 10 7 K, so that 
the pressure profile is proportional to the density profile p = c 2 p/7, where c = ( , yp/p) 1 ^ 2 = 
(jhT / finip) 1 ^ 2 is the sound speed, with symbols h aving their usual meaning, and 7 = 1- 



Radiative cooling is included by using table 6 from Sutherland fc Dopital ( 119931 ). We use a 
time invariable spherical gravity field calculated from the hydrostatic equilibrium at t = 



. . c 2 1 dp 3 c 2 r 

g (r) = — = 2 

7 p dr 2 7 t-q 



i + i- 



-i 



(2) 



To quantitatively analyze the thermal evolution of the ICM we marked several regions 
of the initial ambient gas in the (w, z) plane. Each such region is actually a torus due to our 
2.5D grid. These are called the traced regions. The tracing is done by defining artificial flow 
quantities called 'tracers' in the PLUTO code (IMignone et al.l 120071 ) . that are frozen-in to 
the flow. At t = the tracer of the region centered at (zu, z) = (a, b) and having a radius of 
0.25 kpc is set to £ = 1. Mixing of the traced gas with the ICM or the jet's material changes 
the tracer to values below 1. Since this quantity is advected with mass, the summation over 
all zones % of mass multiplied by tracer value M{ x ^ is constant with time. This was verified 
in the post-simulation analysis. Using these tracers we define the average property Q of the 
traced gas that was centered on (w, z) = (a, b) at t = 



Q 



ab 



(3) 



where Q can be the temperature, the pressure, or the position of the traced gas, i.e., its 
center of mass. 



3. THE GENERAL FLOW STRUCTURE 

In figure [1] we show the density (left) and temperature (right) maps at the end of the 
jet's injection phase (t = 20 Myr). In this simulation the combined two jets (here we simulate 
only one jet) power is Pij = 2 x 10 44 erg s -1 , their initial velocity is Vj = 9600 km s -1 , and 
mass loss rate into the two jets is M.2j = 7M yr _1 . Clearly identified in the figure are 
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the forward shock running into the ICM, multiple sound waves (clearly seen in the velocity 
map presented in figure |2J), and the reverse shock, where the jet is shocked. The shocked jet 
material forms a hot low density bubble, of size ~ 10 kpc at t = 20 Myr. The bubble is the 
region of low density (dark blue and light blue) in figured] (left), corresponding to the hot 
regions (oran ge, red and yellow) in fi gure U (right). The bubble obtained here is similar to 



the results of ISternberg et al. I (120071 ). Such bubbles are observed as X-ray deficient bubbles 
(cavities). Also marked is the contact discontinuity separating the shocked jet's material 
from the shocked ICM gas. At later times mixing prevails and the contact discontinuity 
cannot be identified (see next section). The low post-shock velocity of 85 km s -1 in the 
equatorial plane is because the forward shock in the equatorial plane at t = 20 Myr is very 
weak, having a Mach number of M eq (20) = 1.1. 

To demonstrate the similarity to observed X-ray deficient bubbles, we produce a syn- 
thetic X-ray image by integrating the density squared along the line of sight perpendicular 
to the (x, z) plane, where the x axis on the plane of the sky coincides with the w axis of 
figure [TJ A simulated X-ray image of the cluster inner region at t = 20 Myr is presented in 
figure [31 

The thermal content of the bubble at t = 20 Myr is presented as filling factor (per keV) 
as function of temperature in figure HI The bubble here was defined as material having a 
temperature TJ, > 3T = 10.3 keV, where T = 4 x 10 7 K is the initial ICM temperature. 
At t = 20 Myr the bubble volume is equal to that of a sphere of radius ~ 7 kpc. The 
volume changes a little for small changes of the temperature threshold. The post shock 
temperature of the jet (for 7 = 5/3 and Vj = 9600 km s _1 used here) is T p j = 110 keV. The 
temperature of the bubble is much lower, as energy is transferred to the ICM, by performing 
work on the ICM to inflate the bubble and by heating ICM gas that is mixed into the 
bubble (see section IH). The res ulting volume filling fraction is well below the upper limit 



set by ISanders fc Fabian I (120071 ) for the Perseus cluster, although the work presented here 



is general and we have not set out to reproduce a specific observed system. 

We pick several thin tori at t = and follow the evolution of the material inside them. 
These are called the traced regions. The four traced regions marked on figure [1] were located 
initially at {w,z) = (3,1), (4,2), (6,0), and (9,1); these are called TR31, TR42, TR60 
and TR91, respectively. By t = 20 Myr they moved to (w,z) ~ (5,1), (5.5,2.5), (7.5,0), 
and (10, 1), respectively. We note that TR31 and TR42 were pushed outward and lost their 
circular cross section due to interaction with the shocked gas, both from the jet and the 
ICM. TR60 and TR91 can be seen to have moved a more-or-less uniform distance of 1 kpc 
away from the jet axis. Their motion was set by the forward shock. 



Figure [5] shows the evolution of the flow structure. Already at t — 5 Myr a clear bubble 
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Fig. 1. — The flow structure at t = 20 Myr in the (w, z), where z is along the symmetry 
(jet) axis and w is in the equatorial plane. The simulation is 3D with axisymmetry imposed, 
namely, a 2.5D simulation. Some prominent features of the flow are marked. Left: density 
map with the traced regions marked (TR31, TR42, TR60 and TR91). Due to our 2.5D grid, 
each traced region is a thin torus. Density scale is in units of logp(g cm -3 ). Right: the 
temperature map in units of logT( K). 




to [kpc] 

Fig. 2. — The velocity map in the {w, z) plane at t = 20 Myr. The forward shock, reverse 
shock, and particularly the sound waves, are clearly visible. Velocity scale is in units of 
log v (km s -1 ). 
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Fig. 3. — The projected X-ray map in the (x, z) plane at t = 20 Myr, obtained by integrating 
density squared along the line of sight. The x axis on the plane of the sky coincides with the 
w axis of figure [TJ Only one quarter of the plane was simulated, and it was mirrored twice to 
obtain the full image presented here. The bubbles are clearly seen as regions of low-intensity 
emission. Also clearly seen are sound waves. 
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Fig. 4. — The volume filling fraction (per keV) of the bubble as function of temperature. 
The bubble here was defined as material with temperature above 3 times the initial ambient 
temperature. 



-8- 



is seen. At t = 50 Myr the bubble has moved away from the center by inertia and buoyancy. 
Also seen is a trailing region of vortices that drags ICM outward. By t = 80 Myr we can 
still identify the bubble, although it starts to dissipate. The velocity maps at these times 
are presented in figure [6j Most of the bubble's volume is filled with a large vortex. In our 
2.5D simulations the vortex has a shape of a torus, but in a realistic 3D flow we expect 
the presence of a more complex structure of vortices. The evolution of vortices, inside and 
outside the bubbles, play a major role in the heating and mixing as we show in the next 
section. Here we follow the position of five vortices from t = 50 Myr to t = 80 Myr as 
marked on figure Vortices A, B and b merge, vortex C slowly moves outward, while 
vortex D stays at about the same place. At t = 80 Myr the large vortex A, center at 
(zu, z) v = (4, 30) kpc, includes both the shocked jet material (close to the symmetry axis) 
and the ICM. This holds also to the C and D vortices, centered on (zu, z) v = (4, 15.5) kpc and 
(zu, z) v = (1.5, 9) kpc, respectively. iRefaelovich fc Sokerl ( 120121 ) further discuss the evolution 
of vortices and their role in determining the morphology and evolution of bub bles in groups 
and clusters of galaxies (for early studies of vortices in general contexts see iNormanl Il996 
and references therein). 



4. MIXING AND DREDGE UP 

In the previous section the inflation of 'fat bub bles' and the significant role of vortices in 



the evo l ution of bubbles were demon strated (see also lSternberg et al. Il2007t ISternberg fc Soker 



2008a|JbJ; IRefaelovich fc Sokerl l2012h . We now turn to examine the mixing process. 



We first follow the mixing of the jet's material with the ICM and with the traced region 
TR24, i.e., the traced region whose cross section in the (zu, z) plane at t = is a circle of 
radius 0.25 kpc centered on (zu, z) = (2,4), as well as with the TR31 tracer. Figure [7] shows 
the jet's material in colors, the TR24 gas in white contours, and the TR31 gas in black 
contours; contours show where the concentration of each tracer is one percent. The color 
code gives the fraction of the original jet's material at each point. At early times, t < 10 Myr, 
the material near the axis is almost purely of jet origin (light red). Further away from the 
axis, where the shocked jet's gas resides, the so called 'cocoon', gas from the ICM is mixed 
with the shocked jet material, as can be seen by the dark red regions. At t = 20 Myr the 
shocked jet's gas that fills most of the bubble is heavily mixed with the ICM. At t = 20 Myr 
the jets cease, and a region trailing the main (front) bubbles is formed. At t = 50 Myr the 
main bubble is on the upper left, while at t = 80 Myr the main bubble is outside the frame 
of the lower-right panel of figure [7] (see fig. E]). At t = 50 Myr a substantial mixing is seen 
in the trailing volume, as indicated by the yellow and orange colors. This mixing is more 
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Fig. 5. — The density maps in the (zcr, z) plane at four times: t — 5, 20, 50, and 80 Myr. The 
jet is active during the t = — 20 Myr time period. Density scale is in units of logp( g cm -3 ). 
We follow the position of five vortices and mark them on the panels at t = 50 and 80 Myr. 
Note that vortices B and b merge in the time lapse between t = 50 and 80 Myr. 
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Fig. 6. — Like figure [5] but for the velocity and temperature maps. A large and clear vortex, 
centered at (tu, z) v ~ (2.6, 2.6) kpc, is seen in the shocked jet's gas at t = 5 Myr. The 
temperature scale is in units of logT( K), while the velocity vectors are divided into four 
groups by their length, from longest to shortest in km s : 5000 < v, 1000 < v < 5000, 
200 < v < 1000, and 50 < v < 200. 
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significant at t = 80 Myr. Pockets of ICM gas can be seen within these mixed regions at 
these two times. White contours that mark the position of the TR24 traced region clearly 
demonstrate the mixing and dredge-up of the ICM. Mixing seems to a be the major process 
of heating the ICM along the jet's propagation direction. 

We now turn to the mixing of gas perpendicular to the jets' axis. At early times 
(here t — 5 Myr) the traced region TR31 has been pushed by the forward shock and the 
pressure of the shocked jet's material, while still staying intact. At later times the TR31 
initial morphology (cross section in the w, z plane) is violently disrupted by the turbulence 
(vortices) of the shocked jet and ICM gas. By t = 80 Myr the original TR31 gas is heavily 
mixed with the hot shocked gas, and its temperature increases as we show in the next section. 

Another effect of the jets on the ICM is the displacement of the ICM. We follow the 
center of mass of a tracer by taking the quantity Qi in equation ([3]) to be the location of 
the material Qi — » r\. In figure [H] we present the distance r from the center of the centers 
of mass of four tracers. All tracers are pushed outward at early time as the forward shock 
passes through them. Later they can fall back intact, as TR91 does, or be completely mixed, 
a process that can cause the center of mass distance from the center to increase or decrease 
several times. 

The mixing of material from different regions driven by vortices may affect the metallic- 
ity gradient in the ICM. In order to properly examine this, the simulations need to include 
realistic star formation and metal enrichment due to it, and more important, to include 
metallicity-dependent cooling. The reason for the latter is that higher metallicity regions 
will in general cool faster and fall back to the center. This is beyond the scope of the present 
work, and we leave it for a future work. 

5. HEATING MATERIAL NEAR THE EQUATORIAL PLANE 

To better understand the heating processes of the ICM gas residing near the equatorial 
plane we follow the variation of some thermodynamic quantities of several traced regions. 
In figure M we present the relative changes of the average pressure, temperature, location, 
entropy, and the mixing degree of TR31; TR31 is a traced region whose cross section in the 
(w, z) at t = is a circle of radius 0.25 kpc centered on (zu, z) = (3, 1). The average pressure, 
temperature, and location are calculated by equation and AQ = Q a b(t) — Qab(t = 0). 
The entropy (per particle) is calculated from the Sackur- Tetrode equation, using the average 
temperature and average density of the tracer. Also shown is the degree of mixing, which is 
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Fig. 7. — The concentration of jet material £j Ct (fraction of jet's material at each point) is 
shown by the color coded map in logarithmic scale at four times. It is clearly seen that 
the jet's material mixes with the ICM, and as shown in the next section, heats it. We also 
follow the material of two tracers, TR31 and TR24. Here TR31 is the gas that started inside 
the circular region centered on (w, z) = (3, 1) kpc and having a radius of 0.25 kpc, with a 
similar definition for TR24. The black (white) contours show where the concentration of 
TR31 (TR24) is one percent. Most of the TR24 gas is carried with the jet's gas along the 
symmetry axis. The TR31 gas suffers vigorous mixing near the center, and is heated up by 
this mixing. 
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defined using equation Q, so that 

{ s (4) 

As can be seen in figure [HI initially £ = 1, which means there is no mixing. When the 
shock passes through the region, there is a slight decrease in £. This is a numerical effect - as 
the region shifts its location, the zones on the circumference appear to be mixed, although 
the region stays intact. At later times the value of £ further decreases due to mixing. 

The main conclusion from figure M is that the heating is due to mixing. The forward 
shock that runs through the ICM has no significant lasting influence. The forward shock 
hits TR31 at t = 1 Myr, compresses it and heats it. However, the gas re-expands, and at 
t ~ 10 Myr its temperature returns to its original v alue. Its entropy is somewhat higher than 



its original value, but not by much (see eq. 3.7 in lSoker et al.ll200ll ). A significant heating 
of the TR31 gas starts at t ~ 15 Myr when the degree of mixing with hot gas increases (as 
seen by the decreasing value of £). The mixing is with hot shocked jets' and ICM gas; the 
ICM and the jet materials are already mixed as can be seen in figure At t > 50 Myr the 
temperature drops because mixing with cooler ICM medium starts to be more important 
than mixing with hot shocked jet's gas. In figure [TU] the evolutions of the entropy (left) and 
temperature (right) of TR31 with jet and radiative cooling included (the standard run) are 
compared to cases where either the jet or radiative cooling are not included. Because of 
the stochastic nature of the mixing process, small differences in the initial conditions can 
lead to local differences in the thermal evolution. The differences at early times between the 
standard simulation and the one with no radiative cooling is due to this stochastic behavior. 
At late times the tracer TR31 in the standard simulation is cooler than the case with no 
radiative cooling. 

Figure [TT] is similar to figure fTUl but for tracer TR91 and only for the temperature. This 
tracer stays intact till the end of the simulation (t = 150 Myr), although it loses its circular 
cross section. It suffers a small degree of mixing, which accounts for it being somewhat hotter 
at late times than the case where there is no jet. The rapid temperature rises of TR91 are 
due to the forward shock and sound waves that cross the TR91 material. After the passage 
of each sound wave the TR91 gas cools back. The average temperature of the TR91 gas is 
below its initial temperature, but somewhat above the case without jet activity. This shows 
that the efficiency of heating decreases further away from the center. The heating of this 
region will either take place if the jets of the next episode will have a different direction, 
such that TR91 will be closer to the jets' axis, or if the next activity episode will occur 
after ~ 100 — 300 Myr. In this latter possibility the material of TR91 will be closer to the 
center (see figure |SJ), and will be more likely to be mixed with the hot bubble material. At 
t = 150 Myr the center of mass of TR91 is at a distance of 7.5 kpc from the center (fig. [5}, 
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Fig. 8. — Distance from center of the centers of mass of four traced regions as function of time. 
As before, TR31 stands for the tracer that started from the location (tu, z) = (3, 1) kpc, 
etc. The traced regions start to move outward when the forward shock hits them. At later 
times some of the tracers lose their shapes, and mixing over an extended volume determines 
the location of the center of mass of each tracer. 
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Fig. 9. — Variation of the relative change in average pressure, temperature, location, en- 
tropy, and the mixing degree as defined in equation (j2J), of the traced region TR31. Both 
temperature and entropy clearly show that the major heating process is mixing, indicated 
by decreasing value of £, and not the forward shock that hits the boundary of the tracer at 
t — 1 Myr. The decrease in temperature and entropy at late times is due to mixing with 
cooler ICM regions. The cooling of the tracer comes with heating more ICM. 
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and falling slowly inward at a velocity of ~ 0.01 kpc Myr -1 . At this rate, by t ~ 300 Myr 
TR91 will be close to where TR60 was initially, close enough to the center to be heated by 
mixing (see below). In any case, in the cold feedback mechanism some of the gas does cool, 
falls inward and feeds the AGN. So there is no need for 100% heating efficiency 

Figures [T^] and [T3] depict the time evolution of the average temperature and entropy of 
several traced regions. All regions undergo shock heating early on, with diminishing shock 
heating as the region is further away from the center and at a larger angle to the jets' 
direction. However, shortly afterwards the regions cool down due to adiabatic expansion. 
Subsequent heating is not by shocks, but by mixing of the traced region material with hot 
material from the shocked jet and ICM. Mixing is due to the vortices and turbulent nature 
of the jet-ICM interaction process. 



SUMMARY 



We used the PLUTO hydrodynamic code (IMignone et al.ll2007l ) to conduct 2.5D hydro- 
dynamic simulations, i.e., the flow is 3D but with an imposed azimuthal symmetry around 
the 6 = (z) axis, to study the heating of gas perpendicular to the jets' axis. Only one side 
of the equatorial plane was simulated, and the jet was active for a time period of 20 Myr. 
We marked several intra-cluster medium (ICM) regions, the traced regions, and followed 
the evolution of their averaged thermodynamic variables. Each traced regions has an initial 
circular cross section with radius of 0.25 kpc in the meridional plane (w, z), i.e., it is a torus 
in 3D. At late times the traced region are vigorously mixed with the ICM and lose their 
intact structure. 

We reproduced (fig . [TJ the formation of a fat bubble by a slow massive wide (SMW) jet 
( [Sternberg et al. 1120071). and the f ormation of multiple sound waves with a single jet activity 
episode ( jSternberg fc Sokerl 120091 ). Vortices play major roles in the structure of evolution 
bubbles, as is evident also from figures |5] and EJ In the present paper we found that vortices 
play a major role in setting a complex flow structure that leads to a very efficient mixing of 
the high entropy shocked jet gas with the ICM. As evident from figure [7] the mixing is very 
efficient both along the jet's expansion cone (white contours) and near the equatorial plane 
(black contours). The ICM gas along the jet's cone is dredge up to large distances, while the 
material near the equatorial plane spreads and mixes in a large volume in the inner region. 

Our main results are summarized in figures M and [TTJ These are that (i) the heating 
is very efficient near the equatorial plane up to a distance of ~ 10 kpc, (ii) the heating by 
the forward shock wave is very small, and (Hi) that the main heating is caused by mixing; 
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Fig. 10. — Average entropy (left) and temperature (right) of TR31 for three cases as marked 
in the inset: the standard case where both the jet and radiative cooling are included, an 
active jet included but radiative cooling is not included, and a case with radiative cooling 
but no jet. The initial cooling time of the TR31 gas is r coo i = (5/2)nkT/n e n p A ~ 6 x 10 8 yr. 
The main differences between the standard run and that without radiative cooling is in the 
mixing process, where in the latter case the ICM is somewhat hotter at late times. 
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Fig. 11. — Like figure ITUl for the temperature, but for tracer TR91. 
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Fig. 13. — The evolution of the specific entropy of four tracers as marked in the inset. The 
entropy of each tracer is calculated from its average density and temperature. 



-18- 



mixing is depicted by decreasing value of £. These are further demonstrated in figure [TUj 
where the standard run is compared with similar simulations but without either radiative 
cooling or without a jet. The mixing and heating near the equatorial plane continues to 
counter radiative cooling for times of > 10 8 yr after the jets have ceased to exist. 

Our results further have the following implications. (1) The complex mixing will en- 
tangle magnetic field lines and will suppress any global heat conduction in the ICM near 
the center. Namely, the presence of bubbles in cooling flows implies that heating the in- 
ner region by heat conduction is prohibited. (2) The same entanglement process mixes the 
magnetic fields of the ICM and the shocked jets' material. This will lead to reconnection 
of the magnetic field lines, hence will allow for local heat conduction between the mixed 
ICM and jets' gas. We emphasize the necessity to distinguish between the efficient process 
of local heat conduction (scales of < 0.1 kpc) as opposed to t he inefficient global (scales of 



> 1 kpc) heat conduction process (see review by ISokerll2010bl ). (3) The most crucial region 
for heating in cooling flows is the inner ~ 10 kpc. Our results imply that heating by jets that 
inflate bubble is very efficient in this inner region in all directions, along and perpendicular 
to the jets' axis. (4) The vigorous mixing implies that the region within few x 10 kpc is 
multi-phase. Some small regions will be the coolest ones. If they are not heated by another 
jet activity episode, they will cool, flow inward, and feed the AGN. O ur results therefore 



supply further support to the cold feedback mechanism as suggested by iPizzolato fc Soker 



gas and bv more ( 


etailed studies (e.g. 
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2008; 
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2012; 


Sharma et al. 


201 


2: 


Far age et al. 
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The same feedback mechanism that works in group and cluster cooling flows can be 
the feedback mechanism during galaxy formation if the ISM was hot, i.e., at about the 
virial te mperature the re. Namely, a cooling flow might have existed during galaxy formation 
periods (jSokerl l2010al ). We note though that at the time of galaxy formation the volume 
outer to the inner ~ 10 kpc contains much less mass and its pressure is much lower than the 
values simulated here. This implies the following. All traced regions here are accelerated by 
the forward shock to velocities of ~ 100 — 300 km s _1 , and can be further pushed out by the 
shocked jets' gas. However, they slow down because of the interaction with the ICM further 
out, as evident from figure [HI During galaxy formation, on the other hand, this interaction 
can expel a large fraction of the gas outward. The study of this type of interaction during 
galaxy formation is a subject of a future paper. 
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